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43 Abstract 

In this paper, the partial-wave expansion method is applied to describe 
the difference-frequency pressure generated in a nonlinear scattering of two 
acoustic waves with an arbitrary wavefront by means of a rigid sphere. 
Particularly, the difference-frequency generation is analyzed in the nonlin- 



43 



ear scattering with a spherical scatterer involving two intersecting plane 
waves in the following configurations: collinear, crossing at right angles, and 
counter-propagating. For the sake simplicity, the plane waves are assumed 
to be spatially located in a spherical region which diameter is smaller than 
the difference-frequency wavelength. Such arrangements can be experimen- 
tally accomplished in vibro-acoustography and nonlinear acoustic tomogra- 
phy techniques. It turns out to be that when the sphere radius is of the order 
of the primary wavelengths, and the downshift ratio (i.e. the ratio between 
the fundamental frequency and the difference-frequency) is larger than five, 
difference-frequency generation is mostly due to a nonlinear interaction be- 
tween the primary scattered waves. The exception to this is the collinear 
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scattering for which the nonlinear interaction of the primary incident waves 
is also relevant. In addition, the difference- frequency scattered pressure in all 
scattering configurations decays as r _1 lnr and 1/r, whereas r is the radial 
distance from the scatterer to the observation point. 
Keywords: Difference-frequency Generation, Scattering of Sound by 
Sound, Partial-wave Expansion 



i 1. Introduction 



An outstanding feature of the nonlinear interaction of two or more acous- 
tic waves is a generation of secondary waves having different frequencies, 
namely harmonics, sum- and difference- frequency In the presence of an 
inclusion, this generation is enhanced by two physical effects. First of all, 
the incident waves produce a radiation force through nonlinear interactions 
with the ta c luS10n 6 fl Q 1 As . tesu.t, the inc, US ion b set in m ot ion 
emitting waves which frequencies correspond to the components present in 
the dynamic radiation force. In addition, the primary waves (incident and 
scattered) related to the fundamental frequencies interact yielding secondary 
waves. This process is also known as scattering of sound-by-sound in which 



12 sum- and difference-frequency waves are generated 



7 



l9|. 



Difference- frequency generation is present in several applications like para- 
metric array sonar h|, audio spotlight 11], characterization of liquid- vapor 



15 phase-transition 



12], and acoustical imaging methods such as non 



i6 rameter tomography [1 



14 
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16] and vibro-acoustography 



17 



mear pa 



19]. 



Moreover, parametric arrays have been used to produce low-frequency waves 
in wideband scattering experiments [20j. In this case, the scatterer is placed 
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19 outside the interaction region of the incident waves and the scattering is 

20 treated through the linear scattering theory. This is similar to calibrat- 

21 ing parametric sonars based on measurements of the linear scattering cross- 

22 section 2 



23 Investigations of difference- and sum-frequency generation concerning to 

24 spherical and cylindrical scattered waves were firstly performed by Dean- 



III 



22| | . Scattering consisting of nonlinear interaction of a plane wave with a 



26 radially vibrating rigid cylinder 23[ and sphere 24| have also been analyzed. 

27 Moreover, difference-frequency generation in scatte ring of two collinear plane 



28 waves by means of a sphere was previously studied [25( . However, the results 

29 obtained in this study show that the difference-frequency scattered pressure 

30 has singularities in the polar angle of spherical coordinates (i.e. the angle 

31 formed by the position vector and the z-axis). Furthermore, the difference- 

32 frequency scattered pressure only depends on the monopole terms of the 

33 primary waves. Giving this physical picture, a broader discussion is required 

34 on how to handle the singularities and why the information from higher-order 

35 multipole terms of the primary waves were discarded. 

36 Applications of difference-frequency generation in acoustics generally em- 

37 ploy incident beams which deviate from collinear plane waves. This has 

38 stimulated the investigation of nonlinear scattering of two acoustic waves 

39 with an arbitrary wavefront. Our analysis stems from the Westervelt wave 



40 equation 



261 ] . This equation is solved through the method of successive ap- 



proximations in addition to the Green's function technique. Furthermore, 
appropriate boundary conditions are established to garantee a unique solu- 
tion of the Westervelt equation. The difference-frequency scattered pressure 
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44 is obtained as a partial-wave expansion which depends on beam-shape and 

45 scattering coefficients. Each of these coefficients is related, respectively, to 



46 a complex amplitude o 

47 and scattered waves 27 



a partial-wave that composes the primary incident 
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48 The method proposed here is applied to the nonlinear scattering of two 

49 intersecting plane waves by a rigid sphere. The difference-frequency scat- 

50 tered pressure is obtained in the farfield in three incident wave configura- 

51 tions: collinear, perpendicular, and counter-propagating. In this analysis, the 

52 downshift ratio is larger than five. It is worthy to mention that the collinear 

53 configuration of incident waves has been implemented in vibro-acoustography 

54 experiments [l^], while the perpendicular and counter-propagating: arrange- 

55 ments have been experimentally studied in Refs. Q, [lj]] , respectively. To 

56 reduce the mathematical complexity of the model, the incident waves are 

57 assumed to be spatially located in a spherical region. Even though this ap- 

58 proach is not entirely realistic, experimental accomplishment of scattering of 

59 two located intersecting ultrasound beams was reported in Ref. 2q| ] . 

eo The results show that in the collinear case, the nonlinear interaction in- 

ei volving the primary incident waves (incident-with-incident interaction) and 

62 that of the primary scattered waves (scattered-with-scattered interaction) 

63 are responsible for difference-frequency generation. In the perpendicular and 

64 counter-propagating configurations, difference-frequency generation is mostly 

65 due to the scattered-with-scattered interaction. In addition, the difference- 
ee frequency scattered pressure increases with difference-frequency and varies 
67 with the radial distance r from the scatterer to observation point as r _1 lnr 
es and 1/r. A similar result is found in Ref. 22|, though only monopole sources 
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69 were considered. 



70 2. Physical model 

71 Consider a nonviscous fluid with an ambient density po and an adiabatic 

72 speed of sound cq. The fluid is assumed to have infinite extent. Acoustic 

73 waves in the fluid can be described by the acoustic pressure p as a function 

74 of the position vector r and time t. Absorption effects of a viscous fluid can 

75 be readily included for longitudinal acoustic waves (compressional waves). 

76 In this case, the wavenumber of a single-frequency wave becomes a com- 

77 plex number. However, the account for shear wave propagation, which is 

78 supported in viscous fluids, lies beyond the scope of this study. 

79 2.1. Wave dynamics 

so We are interested on describing how a difference-frequency wave is gen- 
si erated in a nonlinear scattering of two incident acoustic waves by means of 

82 a rigid sphere. The scope of this analysis is limited to acoustic pressures 

83 propagating in the farfield. Up to second-order approximation, the farfield 



84 pressure satisfies the lossless Westervelt wave equation 



30] 



as where (3 = 1 + (1/2)(B/A), with B/A being the thermodynamic nonlinear pa- 

se rameter of the fluid. This equation accounts for wave diffraction and medium 

87 nonlinearity. It is worthy to notice that Eq. ([Q) is valid when cumulative ef- 

ss fects (such as wave distortion) are dominant over nonlinear local effects. This 

89 happens when the propagating wave is far from acoustic sources. When the 



wave is observed near to a scatterer, its pressure should be modified to [31] 



4 V c 2 dt 2 

where <fi is the velocity potential. Note that the approximation p = p holds 
for farfield waves. 

Let us assume that the acoustic pressure is given in terms of the Mach 
number e = vq/cq and e C 1 (weak-amplitude waves), where vq is the max- 
imum magnitude of the particle velocity in the medium. Hence, we can 



expand the pressure up to second-order as 



32] 



p = ep {1) + e 2 p {2 \ £<1 (3) 

where p^\ and p^ are, respectively, the linear (primary) and the second- 
order (secondary) pressure fields. In the weak-amplitude approximation (e <C 
1), the primary and the secondary pressures suffice to describe nonlinear 
effects in wave propagation. Now, substituting Eq. (J3J) into Eq. (pQ and 
grouping terms of like powers e and e 2 , one obtains 

v2 - p< " = °' (4) 

V-^),« = -J^. (5) 
cidt 2 ) p 4 dt 2 

These equations form a set of hierarchical linear wave equations. 



(3 



98 2.2. Linear scattering 

99 Assume that two primary acoustic waves of arbitrary wavefront with fre- 

100 quencies U\ and w 2 (w 2 > oji), propagate toward a scatterer suspended in a 

101 host fluid. The total incident pressure due to the waves is given by 

p l = epocKhie-^ + p^e"*" 3 '), (6) 

102 where % is the imaginary unit, p^\ and pi^ are the dimensionless pressure 

103 amplitudes of the incident waves. When the scatterer is placed in the inter- 

104 action region of the incident waves (see Fig. [T]), two primary scattered waves 

105 appear in the medium. Hence, the primary scattered pressure reads 

Ps = epocKp^e-^ + Ps^e-^), (7) 

loe where p Sj i and p s ,2 are the dimensionless pressure amplitudes of the scattered 

107 waves. Therefore, the total primary pressure in the fluid is then p^ = Pi+p s . 

los It is worthy to notice that the quadratic term d 2 p^ 2 /dt 2 in Eq. (jSJ) gives 

109 rise to waves at second-harmonic frequencies 2u\ and 2a;2, sum- frequency 

no u>i + u>2, and difference-frequency u 2 — W\. These frequency components are 

in distinct and do not affect each other. Our analysis is restricted to difference- 

112 frequency component only. 

in By substituting Eqs. (EJ) and (J7J) into Eq. (Tj0), we find that the primary 

n4 pressure amplitudes satisfy the Helmholtz equation 

(V 2 + ^) r in ) =0, n = l,2, (8) 

\Ps,nJ 
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where k n = u n / Cq is the primary wavenumber. 



Pi.ie 




Figure 1: (Color online) Outline of the scattering problem. Two incident waves of arbitrary 
wavefront with amplitudes p^i and Pi t 2 insonify a target. The observation point is denoted 
in spherical coordinates by r(r, 9, ip), where r is the radial distance from the scatterer to 
the observation point, 6 and tp are the polar and the azimuthal angles, respectively. 



lie The incident pressure amplitudes are assumed to be regular (finite) in the 
n7 origin of the coordinate system. Thus, they are given, in spherical coordinates 
us (radial distance r, polar angle 9, azimuthal angle </?) by 



n 



1,2, 



(9) 



119 where J2 lm — > YliLo Y^ m =~v a ni ls the beam-shape coefficient, ji is the spher- 

120 ical Bessel function of /th-order and Y" z m is the spherical harmonic of /th-order 

121 and mth-degree. The beam-shape coefficients can be determined by using 

122 the orthogonality property of the spherical harmonics. Numerical quadra- 

123 ture can be used to compute these coefficients for waves with arbitrary wave- 

124 front [23,122 1 ■ 



125 The scattered pressure amplitudes are given by 



34] 



Ps,n = SnitfHKrWie, ip), n = l,2 



(10) 
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126 where h\ is the first-type spherical Hankel function of /th-order and s™ t is 

127 the scattering coefficient to be determined from acoustic boundary conditions 

128 on the scatterer's surface. 

129 2.3. Difference-frequency generation 

130 The generated difference-frequency pressure is a second-order field in the 

131 Mach number expansion ([3]). Thus, we may express the difference-frequency 

132 pressure as 

p. = e 2 p c 2 p_e-^\ (11) 

where p_ is the dimensionless difference-frequency pressure amplitude and 
w_ = 0J2 — (jJ\. Substituting Eq. fTTTj) into Eq. (jSJ) we find that j3_ satisfies 
the inhomogeneous Helmholtz equation 

(V 2 + k 2 _)p- = fik 2 V, (12) 

133 where fc_ = uj^/cq is the difference- frequency wavenumber and 

V = P*,lPi,2 + P*i,lPs,2 + P* s ,lPi,2 + P* s ,lPs,2, (13) 

134 with the symbol * meaning complex conjugation. The source term V corre- 

135 sponds to all possible interactions between the primary waves which generate 

136 the difference-frequency pressure. 

137 2.4- Boundary conditions 

138 The uniqueness of solutions of Eqs. (jHJ) and (112"]) depend on the acoustic 

139 boundary conditions across the scatterer object boundary. To find these con- 
Ho ditions the physical constraints of the scattering problem should be analyzed. 

9 



141 First of all, the presence of primary and secondary pressures induces the 

142 object itself to move. Consequently, an acoustic emission by the object takes 

143 place in the host fluid, which means further scattering. If both the object 

144 density is large and its compressibility is small compared to those of the host 

145 fluid, the acoustic emission represents only a small correction to the main 
u6 scattering due to the presence of the object in the host fluid 35J. In our 
147 analysis, this correction is neglected and the object is considered immovable, 
us Therefore, the boundary condition for a rigid and immovable sphere of radius 

149 a is that the normal component of the particle velocity should vanish on the 

150 sphere's surface. 

151 The particle velocity given up to second-order approximation is expressed 

152 as 

V = £V (D +5 V 2 ), £ «1, (14) 

153 where and v^ 2 ^ are the linear and the second-order velocity fields, re- 

154 spectively. Thus, for the linear velocity we have • e r | r=a = 0, where e r is 

155 the outward normal unit-vector on the sphere's surface. From the linear mo- 

156 mentum conservation equation po(dv^/dt) = — Vp^\ we find the following 

157 condition for the primary total pressure 



d(Pi,n + Ps,n) 

dr 



0. (15) 



158 This is known as the Neumann boundary condition. After substituting 

159 Eqs. (Q and (flUj) into this equation, one obtains the scattering coefficient 
wo as s 7 ^ = SnicQ, where 

j[{k n a) 

3rd = 7^T -, (16) 

h\ (k n a) 
10 



161 with the prime symbol meaning derivation. 

162 The second-order particle velocity satisfies the conservation equation 



36) 



rhri 2 ) 

Po^ r + V(^) + £)=0, 



(17) 



163 where C = (po/4)D 2 ( - 1 ^ 2 is the Lagrangian density of the wave, with D 2 

164 being the d'Alembertian operator. The function 0^ is the first-order velocity 

165 potential. Projecting Eq. ( fT7l) onto e r at the sphere's surface, one finds 



dr 



lee Now, using the linear relation p^ 



dC 

dr 



Po 



i) 



(18) 



/dt), Eqs. © and (USD, one 



167 obtains the boundary condition for the difference-frequency pressure as 



dr 



dV 



2kik 2 dr 



(19) 



168 2. 5. Green's function approach 

The solution of Eq. ( II 2p can be obtained through the Green's function 
method. Because the normal derivative of the difference-frequency pressure 
is specified on the sphere's surface, the normal derivative of the Green's 
function on this surface should vanish in order to avoid overspecification in 
the method. Thereby, the difference-frequency pressure amplitude is given 
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in terms of the Green's function G(r|r') by [3j 

r , (V 2 - k 2 _)V 



p_(r) = -/3k 2 _ [ V(r')G(r\r')dV + 
Jv 



Ak^k 



2 



k 

2k 



s/.(£L G(r|iW (20) 



leg where S denotes the sphere's surface and V is the volume of the spatial 

170 region from S to infinity. Note that Eqs. (j2J) and (1T91) have been used in the 

171 derivation of Eq. (1201) . The second term in the right-hand side of Eq. (120]) is 

172 related to the second term in the right-hand side of Eq. (j2J). 

173 The contribution of the surface integral for two interacting spherical waves 

174 (monopoles) is found to be fc?_/(fci&2) 2 in Appendix A. In contrast, it will 



175 be shown in Eq. (1251) that the magnitude of the volume integral in Eq. (1201) 
ne is proportional to /3k-/(kik2)- Thus, the ratio of the volume to the surface 
177 integral is k 2 _j{J5k\k2). It is convenient to write the primary angular frequen- 
ts cies in a symmetric way as follows uj\ = Uq — u^/2 and U2 = Uq + u_/2, 
179 where Uq is the mean frequency. Now the ratio between the integrals can 
iso be expressed as f3~ l [{u)a/ujJ) 2 — 1/4]" 1 . Note that uq/uj- is the downshift 

181 ratio. If the contribution from the surface integral is about 0.01 of that from 

182 the volume integral in water, the downshift ratio should be larger than 5. 

183 Therefore, limiting our analysis to downshift ratios larger than 5, we can 

184 neglect the surface integral in Eq. ([20]) . 

185 The volume integral in Eq. (1201) can be split into two regions: a < r' < r 



186 (inner source volume) and r < r' (outer source volume). In Appendix B 

187 the integral corresponding to the outer volume is estimated for two inter- 

188 acting spherical waves. The result shows that this integral is 0(r" 2 ). It 
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189 will be demonstrated that the inner volume integral evaluated in the farfield 

190 k^r ^> 1 is (^(r" 1 ). Hence, keeping only 0(r _1 ) terms in the difference- 

191 frequency scattered pressure, the contribution of the outer volume integral 

192 can be neglected. 

193 The contribution of the second term in the right-hand side of Eq. f[2"Uj) . 

194 i.e. the term related to local effects, in the farfield is 0(r~ 2 ) as long as the 

195 incident waves behave as 0(r _1 ) in the farfield. Therefore, in the farfield the 
we difference-frequency pressure amplitude is given by 



p_(r) ~ -Pkl [ [ V(r')G(r\r'y 2 dr'dn\ (21) 

J a JQ 

197 where dVt' is the infinitesimal solid angle and the integration is performed on 
we the surface of the unit-sphere Q. 

199 In the region r' < r, the Green's function which satisfies the Neumann 

200 boundary condition on the sphere's surface is given by [^j] 

G = ik_ hf\k_r) X i{k-r')Yr{e, <p)Yr{#> <A (22) 



Lm 



where 

X,(fc-rO = 3i(k-r') - j ^- a) h?\k-r'). (23) 

h\ > (k-a) 

After using the large argument approximation of the spherical Hankel func- 
tion 40[ in Eq. ( )22|) . we find the Green's function in the farfield as 

G = — r l Xi(k-r')Yr(6, ^)^ m *(^', </). (24) 



r 

Lm 
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204 Now, substituting this equation into Eq. ( 121 j) along with Eqs. OH]) and (TTU|) . we 

205 obtain the difference-frequency scattered pressure amplitude in the farfield 

206 as 



p_(r, 9, (f) = ^-/_(r, 0, <£>)- 



hk 2 



k_r » 1, 



20? where 



(25) 



(26) 



Lm 



is the difference-frequency scattering form function. The interaction function 
is expressed as 



57' 



E E 



f ( 2Z 1 + l)(2/ 2 + i; 
4vr(2/ + 1) 



W,0 /-»Z,m .mi* »i2 

A O «i,0,i2,0 O ii,mi,i 2 ,m2 U l,il U 2,« 2 

I 4- n (SI) + « , n (IS) 4 S * « , n (SS) 

X 1 Qhl 2 l + S lA@hhl + ^hPhhl + S l,h s 2,l 2 Q hhl j , 



(27) 



where C,'™ i m is the Clebsch-Gordan coefficient, which come from the 



angular integration through the identity 



4l| 



f Y^Y^rrdn = (-1) 



(2fi + l)(2f 2 + l) 
4tt(2/ + 1) 



W,0 f^l—m 
X O 2i,0,22,0°ii,mi,i 



2, m2 



(28) 



The Clebsh-Gordan coefficient satisfies the following conditions 



42] 



m\ + m 2 = m, 
[Z 2 - ill < I < h + h, 



(29) 
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otherwise it values zero. Furthermore, when mi = m 2 = m = 0the/i + 
l 2 + I should be even else the coefficient becomes zero. The cumulative 
radial functions g( ") stands for each possible interaction of the primary waves, 
i.e. incident- with-incident (II), scattered-with-incident (SI), incident- with- 
scattered (IS), and scattered-with-scattered (SS). They are given by 

& = kikjz- f\i{k-r% 1 {k 1 r')j h {k 2 r , y 2 dr\ (30) 

J a 

fffii = kihk- f X i(k-r% 1 (k 1 r')h\l\k 2 r'y 2 dr', (31) 

J a 

g\f 2l = k,k 2 k. f Xi{k-r')hfXkxr% 2 {k 2 r')r' 2 dr\ (32) 

J a 

e g s J = hk 2 k. f X i(k-r')h\ :2 \k 1 r')hil\k 2 ry 2 dr'. (33) 

208 Equations fl25|) and fl26|) along with Eqs. fl30l) - fl33l describe the difference- 

209 frequency generation in the nonlinear scattering of two primary incident 

210 waves with arbitrary wavefront from a spherical target. 

211 In the upcoming analysis, it is useful to decompose the difference-frequency 

212 pressure amplitude following the contribution of each primary interaction as 

213 given in Eq. (j2"7|) . Accordingly, we write 

-(II) . -(IS, SI) . -(SS) / Q ,\ 

V- — V- +P- + P- , (34) 



214 where the super-indexes stand for the interaction of the primary waves and 

215 j)^ IS ' SI - ) = . According to Eq. ffTB"]) each term in Eq. (134")) is related to 

216 the primary pressure as follows: p*iP%, 2 —> pi , (p*iPs,2 + fttiPip) — > pl S ' SI \ 

217 and p* tl p s ,2 ~> P- S) - 

218 We will show later that the scattered-with-scattered interaction provides 

15 



219 the most relevant contribution to difference frequency generation analyzed 

220 here. Thus, let us examine the asymptotic behavior of Qi^iir)^ with /c_r ^> 

221 1. In doing so, we introduce a new variable u — r'/r in Eq. (1331) and then 

@hhi( r ) = kik 2 k-r 3 / xi{k~ru)hf^ (k\ru)h^ {k2ru)u 2 du. (35) 

J a/r 

222 Since the integrand uniformly approaches to the product of the asymptotic 

223 formulas of the spherical functions with large argument in the interval a/r < 

224 u < 1, then in the farfield this integral has can be written 



a/r _ 



sin k_ru 

2 



h\ 1] '(k-a) 



u 



du. (36) 



225 Therefore, 



d+h-h-l 



2 3 [{k_a) 



2 L va/ \h\ 1} '(k-a) 
x [iir - Ei(2ifc_a)] \, k_r > 1. 



(37) 



226 As a result, the contribution provided by the scattered- with-scattered inter- 

227 action to the difference-frequency scattered pressure varies with the radial 

228 distance r as follows 



p- 



ln r Ao 



229 where A\ and A 2 are constants to be determined from Eqs. ( f25i) - ( 1271) and 

230 ( f37j) . The r~ x lnr term happens only in regions containing primary energy 

231 (volume sources). Furthermore, it is known as "continuously pumped sound 

232 waves" , while the 1/r term is called "scattered sound wave" [9| . 
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233 2. 6. Difference-frequency scattered power 

234 The power scattered at difference-frequency is given by 



p_( r ) = £ _Po£of_ I Re |^v_} ■ e r dn, (39) 
2 Jn 



235 where 'Re' means the real-part and the amplitude v_ comes from the difference- 

236 frequency particle velocity v_ = e 2 cov_e _ ^ _ *. In the farfield, cumulative 

237 effects are dominant in difference-frequency generation. Thus, referring to 

238 Eq. (fTTj) we find that v_ ~ — (i//c_)Vp_. After using this result along with 

239 Eq. (J25D and (EED into Eq. ([39]), one obtains 



p-m-^^E^mi 2 - («) 



240 According to Eq. (|38|) the difference-frequency scattered pressurevaries 

241 logarithmically with the radial distance r. This result is also known for two 



242 



concentric outgoing spherical waves 22j. Consequently, the scattered power 



243 given in Eq. (1401) will increase without limit as r — > oo, unless some account 

244 is taken to absorption processes of the primary waves. 

245 2.7. Series truncation 

246 To compute Eq. (|26|) . we have to estimate a priori the number of terms 

247 L_ in order to truncate the infinite series. This is done by performing a 

248 truncation of the incident partial- wave expansion given in Eq. Let L\ 

249 and L 2 be the truncation orders corresponding to the series expansions of the 

250 primary waves (incident or scattered) with frequency u\ and 0J2, respectively. 

251 The parameters L\ and L 2 are related, respectively, to the indexes l\ and I2 
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252 in Eq. ( 12 7p . To determined Li and L 2 , we employ the following rule 43l. |44 1 



L n ~ k n x + c(k n x) 1/3 , n = l,2, (41) 

253 where c is a positive constant related to the truncation numerical precision, 

254 and x is a characteristic dimension involved in the wave propagation. For 

255 instance, x can be the scatterer radius or a linear dimension of an interaction 

256 region of the incident waves. Once L\ and L 2 are established, the truncation 

257 order L_ of Eq. (|2~5]1 is given through Eq. f ]2"§]) as L_ = Li + L 2 . 

258 3. Results and discussion 

259 To illustrate the solution obtained for the difference-frequency scattered 

260 pressure given in Eq. (123]) . we consider a spherical scatterer suspended in 

261 water, for which cq = 1500 m/s, po = 1000 kg/m 3 , and = 3.5 (at room 

262 temperature). The sphere is insonified by two intersecting plane waves which 

263 are confined in a spherical region of radius R. This region is centered on the 

264 scatterer as shown in Fig. [2J The incident wavevectors are denoted by 

265 and k 2 . Yet this model is not entirely realistic, spatially confined plane 

266 waves with fast spatial decay can be experimentally produced by means of 
26? focused transducers 



29] 



The partial wave expansion of each plane wave is given by 39] 



Pi, n = 47r^z'F™*(^,(p n )jK^r)^ m (^^), r < R, (42) 



Lm 



269 where n = 1, 2 and k n is given in terms of (k n , 9 n , <p n ), with 6 n and <p n bein^ 

270 the polar and azimuthal angles, respectively. Comparing Eqs. (]42p and ([H] 
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Figure 2: (Color online) Scattering of two confined plane waves from a sphere of radius 
a. The wavevectors are denoted by ki and k2. The incident plane waves propagate only 
within the spherical region of radius R. 



271 we find that the beam-shape coefficient is given by 

aZ = ^i l Yr{O n ^ n ). (43) 

272 For radial distances larger than R the incident pressure amplitude vanishes, 

273 i.e. p i>n = 0. Hence, the integration interval of Eqs. f[3"0"j) - f[3"2"j) should be 

274 a < r' < R. 

275 The scattering problem can be further simplified by assuming that one 

276 plane wave propagates along the z-axis, thus, ki = k\e z , with e z is the 

277 Cartesian unit- vector along the z-axis. Whereas the other wave travels along 

278 the direction determined by k.2 = sin(6*2)e z + k% cos(6 , 2)e 2 , where e x is the 

279 Cartesian unit-vector along the x-axis. Unless specified, the scatterer ra- 

280 dius is a = 1 mm, the the cspherical region of the plane waves has radius 

281 R = 2.4 mm, the radial observation distance is r = 0.1m, and the mean- 

282 and the difference-frequency are ujq/2h = 1.5 MHz and U-/2ir = 100 kHz, 

283 respectively. Thus, the downshift ratio is fifteen. These parameters are in 

284 the same range as those used in some nonlinear acoustical imaging meth- 
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291 ] . The size factors involved in the scattering problem are 

2 86 k-R = 1, kiR = 14.5, k 2 R = 15.5, A;_a = 0.41, k Y a = 6.07, & 2 a = 6.49, 

287 k-T = 41.88, k\T = 607.37, and fc 2 r = 649.26. The truncation orders are de- 

288 termined by setting the parameter c = 4 (four precision digits) in Eq. (]41|) . 

289 Hence, the truncation orders for p^P, j3^ s \ p- l \ and p^ ss - ) are respectively 

290 given by (L_, L 1} L 2 ) = (69, 33, 36), (48, 33, 15), (51, 15, 36), (29, 14, 15). 

291 The integrals in Eqs. (!30|) - (l33|) can be solved analytically for arbitrary 

292 combinations of the indexes l,l±, and Z 2 . Nevertheless, the number of terms 

293 in the solution grows combinatorially with the indexes. In the present exam- 

294 pie, the analytic solution of the integrals seems not to be practical. Hence, 



295 the integra 

296 method 



45] 



s are solved numerically by using the Gauss-Kronrod quadrature 



297 The directive pattern in the xz-plane of the difference-frequency scattered 

298 pressure given in Eq. (125]) and produced by two collinear plane waves (# 2 = 0) 

299 are shown in Fig. [3J The dimensionless pressures p^ , p^ s ' SI ^ } and p ( L S ' 1 are 

300 also exhibited. The magnitudes of these functions are normalized to the 

301 maximum value of \p\ which is 0.0807. The contribution from p_ ' is small 

302 compared to other dimensionless pressures. In the region 30° < 9 < 330°, 

303 the difference-frequency scattered pressure is dominated by p^. Both p^ 

304 and p_ give a prominent contribution to the difference-frequency scattered 

305 pressure when 9 < 30° and 9 > 330°. In this case, the contribution of 

306 p^ corresponds to 25% of the scattered difference-frequency pressure. The 

307 magnitude of this pressure mostly occurs in the forward scattering direction 

308 {9 = 0°). We notice that as the radius R of the spherical region increases, 

309 the role of p^ overcomes the contribution of the scattered- with-scattered in- 
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310 teraction. The spatial behavior of p^ ss - ) resembles that of the linear scattered 

311 pressure by the sphere as shown in Fig 0b. 

312 The dimensionless pressure p_ is related to a parametric array whose 

313 primary waves are confined in the spherical region of radius R. We can obtain 
3u an approximate solution of the parametric array pressure in the farfield, when 

315 r = re 2 (forward scattering direction 6 = 0°). To calculate the parametric 

316 array pressure we consider the source term in Eq. (1211) as V = e lk - r ' cose ' . 

317 Moreover, we approximate the Green's function in the farfield to 

ik_(r-r' cos g) 

Airr ' 1 ' 



318 Thus, substituting the source term and the Green's function into Eq. (jzTj) , 

319 we find that the dimensionless parametric array pressure is given by 



ik-r 



p^^-f3k 2 R i - , = 0°. (45) 

3r 

320 Using the physical parameters of Fig. [HI we find good agreement between 

321 the this pressure and pi , with relative error smaller than 9%. This error 

322 might be caused among other things by the presence of the scatterer in the 

323 spherical confining region, which is not accounted by Eq. (145]) . 

324 The directive pattern in the xz-plane of the difference- frequency scattered 

325 pressure produced by two intersecting plane waves at a right angle (62 = 90°) 

326 is displayed in Fig. HI The component p^P corresponds to less 1% of the total 

327 pressure and it cannot be seen in this figure. This result is in agreement 

328 with early studies which state that two intersecting plane waves at right 

329 angle do not produce difference-frequency pressure outside the intersecting 
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(a) 



(b) 



Figure 3: (Color online) The directive pattern in the xz-plane of (a) the difference- 
frequency scattered pressure (normalized to maximum value of \p-\ which is 0.0823) gen- 
erated by two collinear plane waves, and (b) the linear scattered pressures. The physical 
parameters used here are r = 0.1m, R = 2.4 mm, a = 1mm, ujq/2tt = 1.5 MHz, and 
w_/27r = 100 kHz. The corresponding size factors are fc_i? = 1, k\R = 14.5, k 2 R = 15.5, 
k_a = 0.41, kia = 6.07, k 2 a = 6.49, fc_r = 41.88, k x r = 607.37, and k 2 r = 649.26. The 
arrows indicate the direction of the incident wavevectors. 

330 region 7]. The term p^ s,SI - ) does not contribute significantly to difference- 

331 frequency scattered pressure. Thus, p^ is responsible for this pressure. 

332 The two mainlobes of the difference-frequency scattered pressure lies on the 

333 forward scattering directions (6 = 0°,90°) of each incident wave as depicted 

334 in Fig.Hlb. Furthermore, these lobes follow the pattern of the linear scattered 

335 mainlobes as shown in|Hb. 

336 In Fig. [51 we show the directive pattern in the xz-plane of the difference- 

337 frequency scattered pressure generated in the scattering of two counter- 

338 propagating plane waves (62 = 180°). The contributions of andj5^ S,SI ^ are 



339 small compared to that from p 



(SS) 



It is known that the counter-propagating 



340 waves weakly interact nonlinearly 50j. Thus, the difference- frequency pres 



341 sure is practically due to p^ S \ The pressure is not symmetric due to a 
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90 1 90 1.5 




(a) (b) 

Figure 4: The directive pattern in the xz-plane of (a) the difference-frequency pressure 
(normalized to 0.0134) produced by two perpendicular plane waves, and (b) the linear 
scattered pressures. The physical parameters used in the evaluation here are the same as 
those described in Fig. [3l The arrows point to the direction of the incident wavevectors. 

342 difference in the incident wave frequencies. The difference-frequency pres- 

343 sure follows the behavior of the linear scattered pressures shown in Fig. [5jb. 

344 

345 The dependence of the difference-frequency scattered pressure with the 

346 radial distance r is exhibited in Fig. [61 The pressure is calculated in the 

347 forward scattering direction 6 = 0°. In all cases, the main contribution to this 

348 pressure comes from p^ s K analyzed here. Note that according to Eq. (I38p . 

349 the difference-frequency scattered pressure varies as A\r~ l lnr + A 2 r~ 1 . 

350 The scattered pressure varying with difference-frequency is shown in Fig. [71 

351 The pressure is evaluated at r = 0.5 m in the forward scattering direc- 

352 tion 6 = 0°. In all configurations, the scattered pressure increases with 

353 difference-frequency. The difference-frequency scattered pressures due to the 

354 perpendicular and counter-propagating incident plane waves have very close 

355 magnitudes. According to Eqs. (f2"5"l) and (13"5j) . the scattered pressure varies 
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Figure 5: The directive pattern in the irz-plane of (a) the difference-frequency scattered 
pressure (normalized to 0.0178) due to two counter-propagating plane waves, and (b) 
the linear scattered pressures. The physical parameters used here are the same as those 
described in Fig. [3] The arrows point to the direction of the incident wavevectors. 



with difference-frequency as u^f(u^), where / is a function determined in 
these equations. Moreover, by referring to Eq. (125]) one can show that the 
difference- frequency scattered pressure diverges when w_ — > 2uo and u\ — > 0. 
Physically the scattered pressure does not diverge, but decays due to atten- 
uation instead. 

It is worthy to relate our analysis with a previous theoretical study on 
difference- frequency generation in acoustic scattering 25(. We have tried to 
draw a direct comparison between this work and the method presented here. 
Unfortunately, we could not reproduce the reference'ss results due to the 
presence of angular singularities in the difference-frequency scattered fields. 
Therefore, no comparison was possible. Furthermore, we did try to explain 



the experimenta 



368 given in Ref. 



48 



results of difference-frequency generation in the scattering 
. In this study, a nonlinear scattering experiment was per- 



369 formed involving two collinear beams and a spherical target. The incident 
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Figure 6: The (dimensionless) difference- frequency pressure magnitude in the forward 
scattering direction 8 — 0° varying with the radial distance r. The physical parameters 
used in here are the same as those described in Fig. [3J 

370 waves are generated by a circular flat transducer. Despite the authors claim 

371 that the incident beams approach to plane waves, the directive patterns of 

372 the linear scattered waves obtained in the experiments do not follow this as- 

373 sumption (see Ref. |49|). Since the scattering does not involve incident plane 

374 waves, a direct comparison of our theory (for plane waves) and the experi 



48 



375 mental results is not reasonable. However, one of the conclusions of Ref 

376 is that the difference-frequency scattered pressure is mostly produced by the 

377 incident-with-incident and the scattered-with-scattered interactions. This 

378 conclusion is also supported by our results. 



379 4. Summary and conclusions 

380 The difference-frequency generation in the scattering of two interacting 

381 acoustic waves with an arbitrary wavefront by a rigid sphere was theoretically 
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Figure 7: The (dimensionlcss) scattered pressure magnitude versus the difference- 
frequency. The physical parameters used in the evaluation here are the same as those 
described in Fig. [3] 

382 analyzed. The difference- frequency scattered pressure in the farfield was ob- 

383 tained as a partial-wave series expansion. The amplitude of each partial-wave 

384 is given by the interaction function S] 71 , which depends on the observation 

385 distance from the scatterer, the beam-shape and scattering coefficients of 

386 the primary waves. The developed method was applied to the scattering 

387 of two intersecting plane waves located within a spherical region. The di- 

388 rective pattern of the difference-frequency scattered pressure was analyzed 

389 in three incident wave configurations: coUinear, perpendicular, and counter- 

390 propagating. In the coUinear arrangement, the incident-with-incident and 

391 scattered-with-scattered interactions provide a more prominent contribution 

392 to the scattered pressure. In all other configurations, the scattered-with- 

393 scattered interaction prevails over the other interactions. The results show 

394 that the scattered pressure increases with difference-frequency Experimental 
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395 evidence of this feature was reported in Ref. [19[. Moreover, the scattered 

396 pressure was shown to vary with the observation distance as r" 1 lnr and l/r. 

397 Sound absorption effects in the fluid were not considered. If only com- 

398 pressional waves are assumed to propagate in a weakly viscous fluid, the 

399 proposed model can readily accommodate absorption effects by changing the 

400 wavenumber k — > k + ia, where a is the absorption coefficient. Attenua- 

401 tion may affect the obtained results here in at least one way. Both incident 

402 and scattered waves at the fundamental frequencies u± and U2 are more at- 

403 tenuated than the difference-frequency scattered wave. Thus, the nonlinear 

404 interaction range of the fundamental waves in a viscous fluid is shorter than 

405 in a nonviscous fluid. Consequently, a less difference-frequency scattered 

406 signal is supposed to be formed in a viscous fluid. 

407 In conclusion, this article presents the difference-frequency generation in 

408 nonlinear acoustic scattering of two incident waves with an arbitrary wave- 

409 front. This study can help unveil important features of acoustic scattering 

410 not dealt with before. 
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4u Appendix A. Surface integral 

415 According to Eq. (1T91 the surface integral in Eq. (1211 is given by 
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416 This integral will be estimated for two interacting spherical waves. Thus, the 

417 source term V is given by 



V(r') 



Ak—r' 



r2 ' 



(A.2) 



From Eq. flMj) the Green's function becomes 



G = k.a (iV ^(0, ¥>)• 



i,m h\ (k-a) 



(A.3) 



419 Substituting Eqs. ( fA~2|) and ( fA~3|) into Eq. fTA~T]) . yields 



2 A?-^ A. 1 



21,2 
2 



'2i + £;_a) 



(r— a) 



Appendix B. Outer volume integral 



The outer volume integral reads 



(A.4) 



loo = 0. 



G oa {r\v')V{r')r n dr'd£t 



where the Green's function is given by 46| 



(B.l) 



Goo = ik- Xi(k-r)hV(k-r')Yr(e, <p)YT(e' t <pT) t (B.2) 



423 with a < r < r'. We assume that the source term is due the interaction of 

424 two spherical waves as given in Eq. ([Oil . By substituting Eqs. flA~2l) and 
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42 5 (|R2|) into Eq. f lB~Tj) . one finds 



(5k. 



■Xo(k-r) 



oo g2ifc_r' 



(B.3) 



After integrating by parts, we obtain 



fciA;. 



■Xo(*-0 



r 



+ 0(r" 



(B.4) 



427 Therefore, evaluating xo{k~r) through the expressions of the spherical func- 

428 tions, we find J™ = 0(r~ 2 ). 
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